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Abstract 

Bigravity is a natural arena where a non-linear theory of massive gravity can be formu- 
lated. If the interaction between the metrics / and g is non-derivative, spherically symmet- 
ric exact solutions can be found. At large distances from the origin, these are generically 
Lorentz-breaking bi-flat solutions (provided that the corresponding vacuum energies are ad- 
justed appropriately). The spectrum of linearized perturbations around such backgrounds 
contains a massless as well as a massive graviton, with two physical polarizations each. There 
are no propagating vectors or scalars, and the theory is ghost free (as happens with certain 
massive gravities with explicit breaking of Lorentz invariance). At the linearized level, cor- 
rections to GR are proportional to the square of the graviton mass, and so there is no vDVZ 
discontinuity. Surprisingly, the solution of linear theory for a static spherically symmet- 
ric source does not agree with the linearization of any of the known exact solutions. The 
latter coincide with the standard Schwarzschild-(A)dS solutions of General Relativity, with 
no corrections at all. Another interesting class of solutions is obtained where / and g are 
proportional to each other. The case of bi-de Sitter solutions is analyzed in some detail. 
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1 Introduction 



In recent years, many proposals to modify gravity at very large, i.e. cosmological, distances have 
been put forward with the hope to find an alternative explanation to the observed acceleration 
of the Universe. Indeed, one can imagine that this acceleration is not due to some, yet unknown 
and dark, component of the Universe, but to some failure of General Relativity at large distances 
(see e.g. pQ). 

To obtain such a modification of gravity, one of the simplest ideas - in fact already present in 
the original motivation of Einstein to introduce a cosmological constant - is obviously to give some 
sort of mass to the graviton. However, it has long been known that theories of massive gravity are 
rather problematic. In the linearized regime, a Lorentz invariant mass term for the graviton must 
have the Fierz-Pauli (FP) form [2J, or else ghosts appear in the spectrum. With the FP mass 
term, light deflection caused by static sources differs from that of General Relativity (GR) even 
in the limit of massless gravitons, a pathology known as the van Dam-Veltam-Zakharov (vDVZ) 
discontinuity [3]. 

Among the recently studied proposals to modify gravity at large distances, Dvali-Gabadadze- 
Porrati gravity ([4] DGP gravity) has been the subject of numerous works in particular because 
of its ability to produce a late time acceleration of the expansion of the Universe even with a 
vanishing cosmological constant [5]. Another theory, sharing with DGP gravity the same tensorial 
structure of the graviton propagator, and (for that reason) also the vDVZ discontinuity, is non 
linear massive gravity. The latter theory can be obtained from some bigravity theory, where one of 
the two metrics is frozen (see [B] and references therein). On the other hand, bigravity theories have 
the ability to produce an accelerated expansion of the Universe in a non standard way [7J , one of 
the two metrics being then regarded as some new type of dark energy component ( it is interesting 
to note that there is an integration constant appearing in the effective cosmological constant in 
these models, which reminds one of unimodular gravity [8]). The vDVZ discontinuity, without a 
cure, would suffice to rule out any theory in which it appears. It was however proposed that this 
discontinuity could disappear non perturbatively, i.e. that it would not be seen in exact solutions 
of the theory, but only shows up as an artifact of the linearization procedure [9J. While this seems 
problematic in non linear massive gravity [ITU [TT] , there is some evidence that it does work in 
DGP gravity [12], [13l [H] . It is then of some importance to better understand exact solutions of 
various types of massive gravities, among which bigravity theories. In this paper we concentrate on 
spherically solutions of the latter theories. Indeed, in bigravity theories, non-derivative interactions 
between two different metrics can give a mass to some of the spin-2 polarizations and a large class 
of non trivial exact solutions are known. This is not the case for more complicated models such as 
DGP gravity, and thus bigravity theories provides an invaluable framework to investigate various 
properties of "massive" gravity and related models of large distance modification of gravity. 

Interestingly there is another motivation to take a closer look at bigravity theories. Indeed, 
it has been shown [T51 [IB] that the vDVZ discontinuity can also be avoided when the mass 
term breaks Lorentz invariance. The different phases of massive gravity with explicit breaking of 
Lorentz invariance have been further analyzed in [16] . Some examples with spontaneous breaking 
(due to additional vector field condensates) have also been considered in Ref. [17] (see also 
[IB]). The construction in [T7] exhibits ghosts and tachyons at low momenta, but it is nevertheless 
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phenomenologically acceptable for certain choices of parameters. The question of Lorentz violation 
consistent with cosmological observations has been recently studied in [19J. Such Lorentz violating 
mass terms also appear in the so-called ghost condensate scenario with non trivial cosmological 
consequences [20]. It seems worthwhile, in this context, to explore alternative scenarios with 
spontaneous breaking of Lorentz invariance. A straightforward example is provided again by 
spherically symmetric solutions of bigravity theories. 

In Ref. [21] we studied the global properties of a wide class of bigravity solutions, and here we 
will develop the theory of linearized perturbations around some of them. A particularly interesting 
example is the solution where both metrics / and g are flat, but with different values of the "speed 
of light" . This is the simplest case where Lorentz invariance is spontaneously broken. It is also 
interesting in its own right, since it corresponds to the far field in a wider class of spherically 
symmetric exact solutions, of the Schwarzschild form. 

The paper is organized as follows. In Section 2 we review the basics of bigravity. We show that 
if the interactions between the two metrics are non-derivative, then there are always a broad class 
of solutions in the Schwarzschild- (A) dS family, with Lorentz-breaking asymptotics. In Section 
3, we study the perturbations to Lorentz breaking bi-flat solutions. In Section 4 we consider 
perturbations around solutions where the two metrics are proportional to each other, concentrating 
in the case of bi-deSitter solutions. Section 5 summarizes our conclusions. 

While this paper was being prepared, a related work [22] appeared where the study of pertur- 
bations to bi-flat Lorentz-breaking backgrounds is outlined. Where we overlap, our findings agree 
with those of Ref. [22] (see also [23] and [21]). 

2 Exact Solutions 

Following [23], we consider the action 

S = J d 4 *^ (^ + ^) + / d ^V^f (^ + L f ^+S int [f,g] (1) 

Here Lf and L g denote generic matter Lagrangians coupled to the metrics / and g respectively, 
and subindices / and g on the Ricci scalar R indicate which metric we use to compute it. For the 
background solutions, we shall restrict attention to the case where there is only a vacuum energy 
term in each matter sector Lf = ~Pf, L g = —p g , where pf and p g are constant. The kinetic terms 
are invariant under independent diffeomorphisms of the metrics / and g, but the interaction term 
is invariant under "diagonal" diffeomorphisms^l, under which both metrics transform. 

There is much freedom in the choice of the interaction term in (OQ). For instance, Ref. [25] 
considered a non-linear generalization of the Fierz-Pauli mass term 

s int = - ( -J a 4 x(- g n-fnr - gmr - sDis^ - j^), (2) 

6 In principle, we might also include derivative interactions between the two metrics compatible with the diagonal 
symmetry, but in general these terms yield a ghost in the vector sector and we will not consider them here. An 
interesting possibility would be to preserve the independent unimodular diffeomorphisms in the kinetic terms, in 
which case the derivative coupling is possible. 
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(3) 



u + v 
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As we shall see in the following Section, linearized perturbations around asymptotically bi-flat 
Lorentz-breaking solutions of this particular theory show a singular behaviour. For that reason, it 
will be instructive to look at generalizations of (J2J). The most general interaction potential which 
preserves the "diagonal" diffeomorphism takes the form [B] 



where r n = tr[.M n ], n : 1, ...,4 correspond to the traces of the first four powers of = f^ a g a u, 
and V is an arbitrary function. 

There is also some arbitrariness in the way one introduces matter fields, since one has two 
different metrics at hand. One can at least couple matter minimally to any one of the two 
metrics, in which case the other metric can be regarded as some kind of exotic new type of 
matter. Those two choices correspond to the two matter Lagrangians L g and Lf, of action ([1]), 
where it is understood that the matter fields entering into L g and Lf are different. This opens 
the possibility to have two types of matter, one which feels the metric g and the other which 
feels the metric /. In fact one can imagine more complicated situations in which matter fields 
would be coupled to some composite metric built out of the two metrics / and g. If one wishes 
to recover the standard equivalence principle, one should obviously ask that standard matter only 
couples to one metric, and a minimal choice is thus, e.g., that all matter fields appear say in Lf 
(respectively L g ), while L g (respectively Lf), will be simply given by a cosmological constant. 
With such a choice, matter moves along geodesies of the metric / (respectively g), and, provided 
the solutions for the metric / are the same as in standard general relativity (which turns out to 
be possible as will be seen below), there would be no deviations from 4D general relativity seen in 
matter motion. However, various question arise, should one wish to consider bigravity theories as 
realistic. Among those, the fact that bigravity theories are likely to suffer from the same type of 
instability (at the non linear level) discovered by Boulware and Deser in non linear massive gravity 
[26] [27] (see however [3 [7J). Our viewpoint is here more to use those theories as a toy to study 
properties of more complicated models, than pushing the idea that bigravity theories are realistic. 
Thus we will not further discuss in the following phenomenological consequences of various types 
of couplings to matter. In the aim of discussing the vDVZ discontinuity, we will only consider 
the possibility to have two types of matter coupled respectively minimally to metrics / and g, as 
indicated in the action (pQ). 

Let us then introduce the exact solutions subject of this study. The general static spherically 
symmetric ansatz for bi-gravity can be written as [28] 



where the metric coefficients are functions of r. Note that in general it is not possible to write 
both metrics in diagonal form in the same coordinate system. 




(4) 



g^dx^dx" = Jdt 2 - Kdr 2 - r 2 (d6 2 + sin 2 9 dej) 2 ) 
f iiv dx l "dx v = Cdt 2 - IDdtdr - Adr 2 - B (d8 2 + sin 2 9 def 2 ) 



(5) 
(6) 
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For the potential (J2J), and for D(r) 7^ 0, it was shown in [28] that the general solution is given 

by 

g^d^dx v = {l-q) dt 2 - (1 - q y 1 dr 2 - r 2 {dd 2 + sin 2 6d<j) 2 ) (7) 



where 



f^dafdx" = -Ul - p)dt 2 - 2Ddtdr - Adr 2 - 7 r 2 (dfl 2 + sin 

P 



7, 



A = j(l-q)- 2 (p + f3-q-f3q), (9) 
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^"UJ (l-q)- 2 (p-q)(p + P-l-Pq). (10) 

Here /3 is an arbitrary constant, 7 = 2/3 and the potentials p and q are functions of r. Somewhat 
surprisingly, these potentials turned out to coincide with those of the standard Schwarzschild- 
(A)dS family. Solutions of the form (J7J1ED are called Type I. 

Substituting the Type I ansatz into the effective energy momentum tensors which are obtained 
by varying Sint, one readily finds that these take the form of cosmological terms: 

f _ -2 6S int _ A f _ -2 6S int _ Kg 

V - ^sfr ~ K f Uv " ^~ V=9 Sg^ ~ s 



where 



K f 4 V2 

Kg 4 V3 



P u {3v + 9/3(1 -v)}, (12) 
/3' v {3u - 9(3(1 +u)}. (13) 



Note that the corresponding cosmological constants for the metrics ([315]) are not determined solely 
by the vacuum energies pf and p g . They also contain a contribution from the interaction term 
in the Lagrangian. This contribution depends not only on the parameters ( and u (recall that 
v = 1/2 — u), but also on the arbitrary integration constant /3, as seen from 



A crucial observation [21] is that even without assuming any specific form for the functions 
p(r) and q(r), the Type I ansatz leads to energy momentum tensors of the form (1X11) . Here, we 
would like to clarify the reason for that, and to show that Type I solutions (as well as some 
generalizations thereof) exist also in the generic case (T4|). Indeed, for arbitrary metrics / and g 



/T£ = -2((g IfY Lv6$ -Y^n{M n y v V<A , 
g^Tl = -2((g/f)- v (vVSS + Y^ n {M n y u V^j , 



(14) 
(15) 
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where we have introduced the notation 



dr ni ■■■dr 1 



n, 



where I is the number of derivatives. Moving to the frame where both metrics are diagonal, the 
matrix A4 = f^ 1 ■ g can be put to the diagonal form with eigenvalues Aj. Two arbitrary metrics 
g^u and which are solutions of the vacuum Einstein equations will be solutions for bi-gravity 
if all the r n are constant and the eigenvalues of the matrix 



Y,<MY V V^\ (16) 



entering (11411151) are all equal to each other. Note that for a given ansatz, the constancy of the 
traces (or of the eigenvalues) is a frame independent notion. The equations of motion will be then 
satisfied provided that 



(17) 



A, = -2K f C(g/fT (vV - \ 5>r n V^j + K fPf , 

A g = -2K f C(g/f)- v [vV + X - J>r n V™j + K gPg , (18) 

where Af and A g are the cosmological constants. 

Remarkably, the non-trivial background ([THE]) has the property that the eigenvalues of Ai are 
constant 

Ai = {7~ 1 ,7" 1 ,7" 1 ,/97" 1 }, 

which implies 

r„ = 7~"(3 + D, det[M]=f3 1 -\ 

Thus, it is enough to impose 

n 

In the frame where M. is diagonal the previous combination is a constant diagonal matrix with 
only two different constant eigenvalues 



(_ n n 



717 



y{n) I _ 



Both eigenvalues will coincide when 

n 7 - n (-l + f3 n ) V {n) = 0. (19) 

n 

This tells us that for any potential there will exist non-trivial solutions with certain 7 and /3 
satisfying ffl9|) (note that the values of V^ 1 ' depend also on (3 and 7). In addition, (I17H18I) must 



7 



be satisfied. These are three equations for the parameters A/, A g , (3 and 7. Therefore, it is clear 
that one of the effective cosmological constants can be chosen arbirarily. It has the status of an 
integration constant. 



Another interesting class of solutions is obtained by taking / and g proportional to each other, 
but otherwise arbitrary 

Uu = l(x)g^- (20) 

In this case, the matrix M. is proportional to the identity A4£ = 7~ 1 <5^ and the energy-momentum 
tensors (I14H15P read 

Afss = n f r a Ti u = -2(K n - 4u (vv - x> 7 - n y {n) ) k 

KK = H9 m Tl = -2CS7 4 ^ (uV + X>^™ F H S »- 



(21) 



Thus, for any matter content this term just adds to the vacuum energy. From Bianchi identities Af 
and A g must be constant, and / and g must then be solutions of the vacuum Einstein's equations. 
Generically, the expressions for Aj jff depend on 7, so that they imply a constant 7. In this case, 
the parameter 7 is determined through Einstein's equations by noting that (1201) implies 

Rg = iRf- (22) 

Clearly, this class will include solutions in the Schwarzschild-(A)dS family, although non-spherically 
symmetric solutions are possible as well. Note also that such solutions can easily be generalized 
to multigravity theories by deconstructing 5D metrics with a warp factor [29] . 

Let us now turn to the study of perturbations to some of these background solutions. 



3 Perturbations around Lorentz-breaking biflat metrics 

In a theory with two metrics with Einstein-Hilbert kinetic terms and no interaction, there are 4+4 
ADM Lagrange multipliers. When we add a non-derivative interaction which preserves diagonal 
diffeomorphisms, 4 combinations of these may in principle appear non-linearly in the action [6]. 
For these, their equation of motion relates them to the other variables, but they do not lead to 
further constraints. Thus, we have a minimum of 4 and a maximum of 8 Lagrange multipliers for 
20 metric components. Hence, we generically expect a maximum of (10 — 4) + (10 — 8) = 6 + 2 = 8 
degrees of freedom and a minimum of (10 — 8)x2 = 2 + 2 = 4. In a Lorentz-invariant context, 
the first possibility corresponds to a massless and a massive graviton, whereas the second would 
correspond to two massless gravitons. In the Lorentz breaking context, it is possible to have a 
massive graviton with just two physical polarizations [30l [3Tj . 

Let us consider a general potential V[{r n }] as in (jlj). As we showed in the last section, the 
vacuum energies pf and p g can be tuned so that the previous potential has asymptotically bi-flat 
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solutions. At large distances from the origin, these take the form 

SV = V^u, f»v = nffjfw, (23) 

where 

= V - ^S°X (24) 

and n^y = diag(l, —1,-1,-1). The parameters 7 and (3 are related by Eq. (|T9|) . For /3 ^ 1 0, 
we cannot simultaneously write both metrics in the canonical form rj^, and Lorentz invariance 
breaks down to spatial rotations. It will be convenient to introduce the general perturbation in 
the form 

r^T^or+v")' ( 25 ) 

g»u = V^ + h 9 ^, (26) 

where ff v is the inverse of f]^. The perturbation to the metric / has been defined with the upper 
indices, just because this simplifies the manipulations which yield the action quadratic in the 
perturbations shown below. For the remainder of this section, all space-time indices will be raised 
and lowered with the canonical Minkowski metric rj^. The interaction Lagrangian quadratic in 
perturbations then reads 

Lint = L int — \/—gp g — V 1 —j Pf = 

-X{ n2 ^ + V'H/'", + V) + n°(*'oo + rV)( A 'oo + ^ V°) 
-2n 4 (n + /T V°)(^ + h«) + n,{h\ l + /*/) 2 }, (27) 
where, after imposing ffT9l) 

M4 = 4C (j) ' 7l o = 3n 3 -2n 4 + n 2 + 7^^n T -(-l+^)K ) (n) 
n 2 = - ^ n 2 7 " n K) (n) , n 3 = «i/V + ^ n[v - u^V^ - ^ nm 7 -(" +m ^ (ri ' m) , 

n m,n 
ni = n + P^ fe ni -n { _ 1+ pn )V w\ (2g) 

For the sake of simplicity, we will restrict to potentials V^[{t:„}] for which Eq. ffT9l) is independent^ 
of j3, and determines 7. From equation (|2"5|) . this implies n = n 4 . In particular, this class 



7 For P — 1, we have proportional flat metrics the perturbations of which can be obtained from the flat space-time 
limit of the calculations done in the next section. 

8 The case where (TIT)]) is satisfied independently of f3 or 7 leads to the condition 

3n 3 — 3n + n 2 = 0, n 4 = n , (29) 

which, as we shall see, corresponds to the case of no corrections to the Newton's law. An example of an interaction 
where these conditions are satisfied is a potential which is only a function of the ratio of determinants of / and g; 
that is V [{t„}] = V[f / g\. In this particular case, there is an enhanced symmetry under independent "non-diagonal" 
unimodular diffcomorphisms, which do not change the value of the determinants of the respective metrics. 
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includes the interaction (j2J), which, as we shall see, leads to a rather pathological behaviour for 
the perturbations. On the other hand, it is general enough to be representative of generic choices 
of potentials. 

In Refs. [T5l [16] the case of a single graviton with a Lorentz violating mass term has been 
discussed. For comparison with those references, it will be useful to introduce 

wig = — cn , m\ = 0, ml = cn 2 , ml = —cn 3 , ml = —cn^, 

where c > is an irrelevant constant which has the dimensions of mass squared. 

Note that the components h 9 Qi and h f 0t are absent from (1271) . As noted in [22] the absence 



of such terms is a consequence of invariance under diagonal diffeomorphisms in this background. 
In the case of a single graviton (with a Fierz-Pauli kinetic term), the absence of hoi in the mass 
term leads to a very interesting behaviour [16 1 l30 | 132]. where the two polarizations of the massless 
graviton acquire mass, while all the other modes do not propag ate@ 

Let us now investigate whether a similar phenomenon occurs in our model. The situation is not 
directly reducible to that of a single graviton, since the equations of motion are not diagonal. Also, 
the kinetic term breaks the Lorentz invariance. It is convenient to decompose the perturbations 
into irreducible representations of the spatial rotations, 

h X oo = 2^4 X , 



h x 0t = B x + V x 



h% = 2^ x 6 tJ - 2E X - 2F ( x j) - t$ , (30) 

where t x u = t x ^ i = V x ii = F x ii — for X = f, g, and all space-time indices are raised and 
lowered with the metric rj^. 

To second order in the perturbations, the kinetic terms in ([1]) can be written in terms of these 
scalar, vector and tensor variables as: 

l k = JL{ - V: \t% - \ (vf + Pf) a (y? + if) + 4A^ {a? e») 

-2r^r - e(^) 2 } + ^{ - - ^— (yi + pM) a (V + 

+4/T 1 A?/> / (A f - (3B f - (3 2 E f ^j - 2tP f AtP f - 6/3(^) 2 } (31) 

where □ = rf v d il d vi kf = r )~ l f3 1 ^ 2 Kf and dot means a derivative with respect to time. At the 
linear level, the transformations generated by independent diffeomorphisms 5x fl = in each one 
of the metrics can be expressed as 



9 It should be stressed that the absence of Oi components is a peculiarity of the background considered. By 
suitable adjustment of the vacuum energies, the theory we are considering also admits the Lorentz preserving 
vacuum of type II, where / M „ = = rj^ v . In that case, the interaction term leads to the Fierz-Pauli mass term 
for a combination of the two gravitons. This mass term does contain the Oi components. 
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Note that the kinetic term is written in terms of the following quantities: 

tf p Vf + Ff, V, .V"- />' -/>• 
tl, V/ + (3Ff, if,, A f - p& - (3 2 Ef, (33) 

which are invariant under both gauge symmetries. On the other hand, the full action (including 
the mass terms), is invariant only under the diagonal gauge symmetry 

% = (34) 

No second order scalar combination of h x 0i is invariant under this gauge symmetry, which implies 
that those terms are always absent (cfr. (1271) ). We may now analyze the propagating degrees of 
freedom. 

3.1 Tensor Modes 

The linearized Lagrangian for the tensor and vector modes can be expressed as 

1 r 1 „„ „ 1 



L ^2^{-^-2^ + F ") A ^ + F ")} 



Upt{ 3 - ^ (yl + pf?) a (V/ + ppf) } 



2k, I 4 

-^{^(fa + - 2n ^ F " + F D A ( F i + F !)\ ( 35 ) 

where Kf = 7~ 1 /3 1//2 K/. The corresponding equations of motion in Fourier space read 

Jt'L k% • K g M±n 2 {t% • t£) (36) 

(3u% = k% + KfM 4 n 2 (tfj + t{ 3 ) (37) 



t{ 3 k% ■ k ; .u'// 2 (/;' ; • //; 

from which we obtain the dispersion relations 



<4 = 4 ((Z 3 + !) k2 + R oM 4 ± v / ((/3 + l)k 2 + /t M 4 ) 2 - 4/3k 2 ( Kl M 4 + k 2 )) (3? 



where k = n 2 (j3K g + kf) and «i = n 2 [K g + Kf). 
At high energies, we have 

^k 2 , c^w/THe 2 . (39) 

In this limit, each one of the two gravitons propagates in its own metric (with the corresponding 
"speed of light" 0) along null directions fc M = (a>, k) satisfying 

9x kyk v ~ 0. 



10 Super luminal propagation has previously been considered in several contexts (see e.g. [33] for a recent discus- 
sion). Clearly, such propagation cannot by itself be considered pathological. Indeed, in the present case we always 
have superluminal propagation from the point of view of one of the metrics, whereas there is not any supcrluminal 
propagation from the point of view of the other metric. Nevertheless, the global structure of non-linear bi-gravity 
solutions is complicated in general, and its interpretation is far from trivial [21, 23 . 
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The low energy expansion of (1381) is given by 



^ = ^ k 2 + (k 4 ), (40) 

-i=^+( ;/ +/? > l k 2 + o(k 4 ). (4i) 



/3 \(3k f + (3 2 K, 

The first dispersion relation corresponds to two massless polarizations which propagate at the 
"intermediate" speed 

,2 



U0 2 _ Ki K g + Kf 



k 2 k (3n g + kf 

Note that for j3 > 1 we have < c 2 < 1, while for (3 < 1 we have 1 < c 2 < The second 
dispersion relation, Eq. ( I4T1) . corresponds to two massive polarizations. It is easy to check that 
the graviton polarizations are stable and tachyon free as long as > 0, in the whole range of 
momenta k. The second dispersion relation (|4T|) corresponds to the massive graviton. 

3.2 Vector Modes 

From the Lagrangian (1351) . we find that Vf and V/ do not appear in the interaction term. Varying 
with respect to the vector fields we have, 

A(Vf + F°) = (42) 
A (yf + if) = -M 4 n 2 K g A (F? + if) (43) 

A(V/ + [3F/) = (44) 
A (v/ + fiPty = -M 4 n 2 k f A (F? + if) . (45) 

We can always use the diagonal diffeomorphism invariance to work in the gauge where Vf = 0. 
It then follows from (1421) that Ff = Fi(x) + ff (t), where Fi are arbitrary functions of position 
and fi are arbitrary functions of time. The latter are in fact irrelevant, because F t x enters the 
metric only through spatial derivatives. Formally, we may describe this as a gauge symmetry 
F t x — > F x + f x (t), which we can use in order to write, without loss of generality, 

Ff = F t {x). 

It then follows from fj4*3l that 

where again we eliminate the additive time dependent part. Finally, from (1441) we obtain 

v/ = m 

where fi are new arbitrary functions of time. This is not a desirable situation, since it means 
that the initial conditions do not determine the future evolution of V/ . Technically, the absence 
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of the fields V 9 and V/ in the mass term leads an enhanced gauge symmetry in the linearized 
Lagrangian. Indeed, we can consider independent gauge transformations for each of the metrics 

of the form £f = (,f(t). As we have discussed, these do not affect the , but can be used to 
give both of the an arbitrary time dependence. 

3.3 Scalar Modes 

The Lagrangian for the scalar modes can be expressed as 

L s = ^{2A^ 9 (A 9 - B 9 - E 9 ^j - i) 9 Ai\) 9 - 3(</> 9 ) 2 } + —{2/3-^^ (A f - (3B f - (3 2 E f ^j 

-ip s Aip s - 3/?(^) 2 } - ^-{n 2 {3(4> g + ^ f f + (A(E 9 + E f )f - 2(^ 9 + *p f )A{E 9 + E f )} 

+n {{A 9 + /TW) {A 9 + p- l A f - 2[3(t/; 9 + ^) - A(E 9 + E f )})} + n 3 {3(^ 9 + *p f ) - A(E 9 + E 1 )} 2 }. 

Let us first study the non-homogeneous modes. The mass terms do not depend on B 9 nor on 
so those fields are Lagrange multipliers, just as in Einstein's gravity. Variation with respect to 
these fields yields 

Aip 9 = Aip f = 0. (47) 
The variation with respect to A 9 and A^ yields the constraints 

A 9 = -p~ 1 A f + 3(ip 9 + ip f ) - A(E 9 + E f ) + — Aip 9 , 

M 4 n K g 

r = ^f + f(t). (48) 
K f 

Once we substitute the first of these constraints in the Lagrangian, the quadratic term in E h and 
E l takes the form 

(n 2 -n + n 3 )(E h + E 1 ) 2 . (49) 

We can now distinguish two different cases, neither of them with propagating scalar degrees of 
freedom. First, if the coefficient n 2 — n^+n^ does not cancel, the equations of motion for E h and E 
result in a new constraint which determines these fields, and upon substitution into the Lagrangian 
we are left without any scalar degrees of freedom. If the coefficient cancels, as happens for the 
potential (0), E 9 and E f are Lagrange multipliers appearing in the gauge invariant combination 
E h + E . After using (I48p . the variation with respect to E yields 

Aip 9 = Aip f = 0. (50) 

The Lagrangian cancels after substitution of these constraints, and there are no propagating 
degrees of freedom. Note that in this last case the combination E h + E l , is not determined by the 
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equations of motion. Again, this is not a desirable feature, since it means that the value of this 
combination, which is gauge invariant under the diagonal diffeomorphisms, is not predicted by 
the linear theory. Nevertheless, we expect that higher order terms in the expansion will determine 
E h + E l , since there is no symmetry in the non-linear Lagrangian under which this quantity can 
be "gauged" to arbitrary spacetime dependence. 

Concerning the homogeneous modes, after using the constraints we are left with two modes -0 
and ip 9 which have a negative definite kinetic term. Nevertheless, the dispersion relations for the 
degrees of freedom which diagonalize the equations of motion are uj 2 = and uj 2 = M 4 n 2 (kf+n g ) > 
0, so there is no classical instability associated to these modes. 

3.4 Coupling to Matter and vDVZ discontinuity 

All known explicit and non-singular exact solutions of bigravity (which we reviewed in Section I) 
are also solutions of GR. This immediately suggests that the vDVZ discontinuity may be absent 
altogether in this theory at the non-linear level. Also, from the analysis of perturbations done in 
the previous Section around the Lorentz breaking background, it is clear that the situation here 
is very different from that of ordinary massive gravity. The massive spin-2 graviton has only two 
physical polarizations (as opposed to the five polarizations of the ordinary FP massive graviton), 
and there are no propagating vectors or scalars. 

Let us consider the coupling of the linearized theory to conserved sources. To this end, we 
introduce the couplings 



(51) 



where Tg u and T? v are conserved, i.e. d^T^ u = and r] plJ ri p0i d a T^ 1 ' = 0. In terms of the decompo- 
sition flSUl). we have 



S ma tt = ^Jd 4 x (••/•///;;' + 2T°\Vf + Ff) + 2T 9 0C V + 2Tfr) 

+ k / d 4 x [-ij/iy + 2Xf{V/ + pP() + 2T™& + 2T)V / ) • (52) 

where we have introduced the gauge invariant combinations 

$ 9 = A 9 - B 9 - E 9 , ¥ = A f - (3B f - (3 2 E f . 

Inverting the equations of motion for the tensor modes in the presence of the source T % \ we find 

9 = A 9 (k 2 - (3u 2 + k f M*n 2 )T* - A /S M 4 n 2 ^ 

ij LU 2 {(3tu 2 -(k f + (3K g )M^n 2 } + k 2 {(k f + K g )M^n 2 -((3 + l)uj 2 } + k^ [ ' 

and an analogous expression for t^: 

f A/ (k 2 - uj 2 + K g M 4 n 2 )Tl 3 - XgkjM^TPj 

tij ~ u^^ 2 -(k f + f3K g )M^n 2 } + k 2 {(k f + K g )M' i n 2 -((3 + l)uj 2 } + k^ { } 
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In the limit M 4 — > this reduces to the standard expression for linearized GR. 
For the vector modes, the equations of motion read 

A(Vf + F?)=\ g K 9 T™ 
A (Vf + F?) = -M 4 n 2S A (if + F/) + X g K g ff (55) 
A(V/ • IFf) Xj iH/lf 
A (y/ + (3Fty = -KfP-WnzA (if + if/) + XfKffXif. (56) 

It follows immediately that A(if + if) = 0, and therefore the term proportional to M 4 vanishes. 
This means that there is no difference with the GR results for each one of the metrics. 

For the scalar part, we may start with variation with respect to Bf, which yields the constraints 

Cx = (57) 

where 

C g = 4A^ 9 + A g /%T°°, C f ee 4A^ / + Xfh f (3Tf. 
Variation with respect to A x gives 

Cf = C g (58) 

and 

C + = C f + C g = 2M\k f + K g )(A + - + AE + )n , (59) 

where A + = A g + f3~ l Af, ip + = ipf + ip g , and E + = Ef + E g . Variation with respect to Ai^ x 
yields, with the help of (i5T|) . 

n A + = (n 2 + 3n 3 )</> + - (n 2 + n 3 )AE + . (60) 
Substituting into (1591 . we have 

C+ = 2M 4 (k f + K g )[(n 2 + 3n 3 - 3n )ip + - (n 2 + n 3 - n )AE + ] 
and using (!57j) . we have 

4(n 2 - n + n 3 )A 2 £ + = -(n 2 + 3n 3 - 3n )(A / /i / /3T / 00 + A fl /c fl t°°). (61) 

For (n 2 — no + n 3 ) 7^ 0, this determines E + in terms of the sources. The solution will depend on 
an arbitrary time independent mode E (x). 

For the singular case (n 2 — n + n 3 ) = 0, Eq. (1611) do not determine E + at all. Instead, it 
imposes some non-trivial equations to be satisfied by the sources, 

XfKfPlf = -\k 9 T™ (62) 

which seem hard to motivate. Thus, coupling to the sources seems rather inconsistent in this case, 
unless (n 2 + 3n 3 — 3no) = as well. But this would imply n 2 = 0, in which case the tensor modes 
are massless. 
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In the generic case, the solution for the ip potentials is of the form 



Aip 9 
Aip f = 



4 ~ 9 

KfXfP Q n 1 



(63) 



where C + (x) is entirely determined by initial conditions. 

Finally, variation with respect to i/jf and ip g leads [after use of (|63|) ] to the following equations 
for the gauge invariant potentials: 



A$ s 



P f /3T°° + T; - ^/3 2 T> 00 J + + K f M A n 2 AE+. 



(64) 
(65) 



where 
AE, 



1 



^2 + ™3 ~ TIq 



n 2 + 3n 3 - 3n 
4A 



(«,A,2f + /C/A//32J - C+) + 



1 



-C, 



2M 4 (kf + Kg, 

(66) 

In general, the solution depends on an arbitrary "initial" function C + (x). This corresponds to a 
mode with dispersion relation u 2 = in the linear theory. It was argued in [16] that in such cases, 



from higher order terms the expected dispersion relation will be of the form ur ~ p , and in this 
sense C + corresponds to a slowly varying "ghost condensate" [34]. In what follows, we shall take 
the initial condition C + (x) = 0. 

For n 2 — Uq + n% ^ 0, the solution is of the form 



A^ 



A 9 ' 



(67) 



and 




A$ 3 = 


K g \g I 




4 V 


A& = 


Rf\ f p 



9 9 /± 9 



K q M^TL 2 



n 2 + 3n 3 - 3n , 00 

\ K 9 A 9 1 g 



+ Tf - —P'T™ ) - 



4A J n 2 + n 3 - n 
Rf(3M A n 2 \ n 2 + 3n 3 



RfXtfTf) , (68) 



4A 



3n 



n 2 + n 3 - n 



{k 9 \ 9 T™ + R f \ f (Hf) 



Hence, there is a well behaved massless limit, with corrections of order M 4 A~ 2 to the gauge 
invariant potentials $ and if). This means, in particular, that there is no vDVZ discontinuity. 
This is quite analogous to the "half massive gravity" model discussed by Gabadadze and Grisa 
[3T] (see also [30J). The additional terms lead to corrections to the Newtonian potential. The sign 
of this correction can be positive or negative, depending on the values of the numerical coefficients 
rti. For isolated sources, such corrections scale like the square of the graviton mass m 2 ~ kM 4 
times the "Schwarzschild" radius r s corresponding to the given source, and grow linearly with the 
distance r. Parametrically, the potential takes the form 



!>jv + m r s r, 
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where 0at is the standard Newtonian potential. Linear theory breaks down at large distances, 
when the second term is of order unity. It would be interesting to try and match this solution to 
a non-perturbative exact solution which is well behaved at infinity. 

As we stated before, the case of no correction to the Newton's law corresponds to the case 
where f lT9|) is independent of (3 or 7 (cfr. ff29l) ). 



Finally, we note that the simple interaction term ([2]) first considered in [2B^ ¥2E\ happens to 
land on the special case 

n 2 -n + n 3 = 0, 

where the above expressions for the gauge invariant potentials are singular. The origin of the 
singularity is the following. After substitution of the constraints fl60|) . the linearized action no 
longer depends on AE + . In particular, the absence of this variable results in the unwanted 
restriction (161Z1) on the sources 1 H Nevertheless, beyond the linear order, the action will depend on 
AE + , and hence the "restriction" will no longer exist. Rather, a nonlinear equation will determine 
the value of AE + . Can we nevertheless try to find classical solutions in a perturbative expansion? 
The above considerations suggest an expansion scheme for the singular case n 2 — + 113 = 0, 
where E + is treated as a much bigger quantity than the rest of the linearized fieldsE^l (such as if). 
Heuristically, the size of AE + can be estimated as follows. Instead of perturbing the flat solution 
Eq. (1231) . we may consider the quadratic action for perturbations around a solution which differs 
from the original by 0(h). The expansion around this new solution will havj^l 

n 2 — n + n 3 = 0(h). 

From (1631) . we have 

Aip ~ kT = A<pN, 

where 0tv stands for the potential corresponding to the given source in Newton's theory. From 
f lBU]) . A$ ~ O(kT) + 0(m 2 AE), where in 2 ~ kM 4 denotes the graviton mass squared. From flBB]) . 
we have AE x 0(h) ~ 0(nT/A). This suggests the hierarchy 

AE > ij), $ ~ max(^, m 2 E). 

Taking 112 — 110 + 113 ~ max($, AE) ~ max(ip, m 2 E, AE) ~ AE (1 + m 2 /A), this leads to the 
estimate 

(AE) 2 1 . 

1 J l+m 2 /A 

1 /2 

For distances shorter than the inverse graviton mass, we have AE ~ > and hence we may 
expect 

$ ~ <j) N + (in 2 / A)^ 2 . (A > in 2 ) 



11 This accidental symmetry is similar to that which exists in ordinary massive gravity where the linear action 
propagates 5 dot whereas a new ghost-like dof appears at the non-linear level [35, 26 . However, in that case the 
accidental symmetry corresponds to a symmetry of the massless theory and no further constraints are needed in 
the sources. 

12 Some of the linearized fields will be of the order of E as is clear from (|60p. 

13 All the coefficients will have corrections of order 0(h). However, for the rest of coefficients one expects that 
they will yield second order small corrections. 
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At distances which are large compared with the inverse graviton mass, AE ~ (Acpj^/m 2 ) 1 ^ 2 , and 
we expect 

$ ~ (m 2 /A) 1/2 0^ /2 . (A < m 2 " 



These very crude arguments seem to indicate that, also in this special case, there is no vDVZ 
discontinuity. However, for finite m, there are significant modifications to the value of the "gauge 
invariant" potential $ which determines the motion of slowly moving particles. For isolated 

1 /2 

sources, such modifications scale like r s , where r s is the "Schwarzschild" radius corresponding 
to the given source. They grow with the distance as r 3 / 2 below the graviton Compton wavelength 
m _1 , and as r 1 / 2 for larger distances. The potential $ becomes of order one for r > wr 2 r~ 1 ^ 
beyond which we enter a non-perturbative regime. It would be interesting to confirm this heuristic 
analysis in a numerical study of a spherically symmetric solution with sources. This is left for 
further research. 

4 Perturbation theory of Proportional de Sitter Metrics 

As stated in Section 2, another interesting class of solutions of bigravity can be constructed from 
two proportional metrics with a constant proportionality factor. Let us define our perturbations 

as 

g^ = ^ + h%,, (69) 
f»v = 1 ~ 1 (W + hf). (70) 

All indices will be handled with the metric. 

We first focus in the interaction term for a general potential (jl]). Using f[2"Tj) we can write 



Lint = a-gn-fyv[{r n }] + + y/-f 

Kg Kf 

= -^V=n{* + h T)( h U + Ku) - m2 s ( h9 + h f) 2 } . ( 71 ) 

where indices are manipulated with the metric f^, e.g. h 9 = Q^h 9 , and 



m 2 = ( -uvV + (u-v)J2 nj~ n V^ n) + ^ nm^ {n+m) V {n ' rn) J , 

\ n n,m / 

m 2 = -A K+ C,i iv n 2 l~ n V (n) . (72) 

n 

We have also introduced an effective Newtons's constant k + for later convenience. 

Note that the massive graviton corresponds to 7t+ = (h g + hf)^. This is to be expected, as 
for h 9 = —hf the metrics are still proportional and therefore the perturbations are standard 
massless gravitons of GR in vacuum. Also, in the present set-up, h^ v are the quantities invariant 
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under the diagonal diffeomorphisms. Notice also that the mass term does not have in general a 
Pauli-Fierz form, 

m*{hl-h^h%). (73) 

This particular form can only be achieved by properly tuning the parameters. This is in contrast 
with other ways of getting massive gravitons, such as dimensional reduction, where the original 
symmetry group is much larger. Here, the degrees of freedom of the original theory are 8 which can 
be split into a massless graviton with 2 polarizations and a massive graviton with 6 polarizations^!. 
The expression of the massless graviton as a linear combination of the metric perturbations will 
be given below. 

From Eq. (ITTj) we note that whenever m t = there is an enhancement of the gauge symmetry, 
which now admits all transformations which leave the traces h g and h f invariant lj. This cor- 
responds to the transverse subgroup of the diffeomorphisms, which has been recently considered 
in [37] . In this special case the gauge symmetry is enough to have just two massless gravitons 



propagating [3 

Let us now consider the case of generic m s and m t . For simplicity we will concentrate on 
perturbations around de Sitter solutions which will be foliated by spatially flat sections, 

tt^dx^dx" = a{r]) 2 {dr] 2 - S^d^dx 2 ), (74) 

where a(r/) = — (Hr))~ l , H 2 = A g /3 being a constant and r\ e (— 00, 0). The kinetic term in ([1]) 
will be given by (cfr. flTTl) ) 

Lk = -^-V=9~ (R g + 2A g ) -J-y/Zf(R f + 2A / ), (75) 

with Af = j~ x A g . To second order in perturbations we can rewrite the kinetic term in terms of a 
massive and a massless field, 



1 



* = -77— (R 9+ + 2A 9 ) - — V^Z (R g _ + 2A g ) + o(h s ), (76) 



where k_ = k + = k 9 k 1 (l + k), with k = ^n g K f l , g-_^ = + and g +lxv = + h+ u 
Besides, we have introduced the massive and massless combinations 



ht» = h% + hl u , h~ u = (1 + k)- 1 - Kh^) . (77) 

The dynamics of the massless part is well known. One easily finds that only the tensor modes are 
dynamical. For the generic massive theory in de Sitter space, studying the longitudinal mode of 



14 The number of degrees of freedom coincides with that of higher derivative gravity |36j . 

15 This happens in the case when the derivative of Eq. p^|) with respect to (3 vanishes at (3 = 1. For the case © 
this amounts to 7 = 2/3. 

16 At first sight, this seems to contradict the results of Ref. [55] . where it is shown that we cannot have two 
massless interacting gravitons. However, the starting point in [38] is a free Lagrangian invariant under linearized 
diffeomorphisms. As shown in [37j . there are Lagrangians invariant under transverse diffeomorphisms which prop- 
agate just massless spin-two particles. An extension of the analysis of |38j to the transverse subgroup is currently 
under investigation [39 (see also |40j). 
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the massive representation we would argue that the only ghost-free possibility is the Fierz-Pauli 
mass term, mf = —m 2 [21 SE]- However, in general, this mode decouples only at high energies 
(larger than a combination of the rest of relevant mass scales). For intermediate energy scales, the 
longitudinal mode is coupled to another scalar mode which can modify this picture [HI [27]. Also, 
the curvature scale H could play a role in making these intermediate scales phenomenologically 
relevant. We will study this possibility directly in the unitary gauge. 

Let us first split the degrees of freedom of the massive combination into scalar, vectorial and 
tensorial modes, 

h+ = 2a(r ] ) 2 A J 

h+ = a{ri) 2 {B >i + V^ 

h± = a(r ] ) 2 (2ijd lj - 2E tij - 2F {hj) - (78) 

where ip, B, and E are the scalar modes, Fi and Vi are vector modes, and £y is a tensor mode. 
The vector modes are divergenceless and the tensor modes are transverse and traceless. 

The expansion of the kinetic term in this foliation can be extracted from the usual expansion 
in de Sitter space (see e.g. |42j, notice however the difference of convention). One finds 

- ajr / d ' x + 2A ^ = - 2ir / d4 ^){^ D ^ + \^ + F i)^ v * + F i) 

+6(V>' + UAf - 2AiP(2A - V) - 4A(5 + E')(ip' + HA)}, (79) 

where TC = a(r])'/a(ri) = a(i])H and the prime refers to derivative with respect to the conformal 
time r). We have also introduced the d'Alembertian □ = rj^d^du and the Laplacian A = didi. 
The interaction term (ITT|) reads 

L in t = ^-a{ V f\m 2 s {A + AE - 3^) 2 

1 / DAD s 

--m 2 t - 2{VM + FiAFi) + A(A 2 + ^— + (AE) 2 + 3^ 2 - 2^AE)jj. (80) 
We can now analyse the different components in turn. 



4.1 Tensor and Vector Modes 

The action for the massive tensorial modes is simply 

{t) 5S 2 = J dx 4 a 2 (r]) (t^Ut^ + a(r/) 2 m 2 ^^) . (81) 

From this equation we can read the mass of the graviton which will be given by m 2 , and the 
tachyon-free condition will simply read 

m 2 > 0. 
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Regarding the vector modes, their action is 

^5S 2 = J dxV(r ? )((V i + ^)A(^ + ^) -a 2 (^)^ 2 (^ + F i AF l )). (82) 

The field V m enters the action without time derivatives, and thus its variation yields the constraint, 

A(Vi + F'A = a^fmlV, = m 2 ( V )V t . (83) 

Taking this constraint into account, the action for the vectorial modes up to second order can be 
written as 

^5S 2 = -L J dVMm 2 W ( F i A _ m{vr F i + W*) (84) 

This Lagrangian has the usual signs, and thus no ghost or tachyons appear in the theory for 
m 2 > 0. More concretely, we can canonically normalize the previous field equation with the field 
redefinition 

F ' =mM ^4A^W) F '- (85) 

We conclude that the only constraint we get form the analysis of the vector and tensor modes is 
m 2 > 0. 



4.2 Scalar Modes 

From (1791) and (ISO]) , the second order Lagrangian for the scalar part reads 
1 



^6S 2 



2ka 



J d 4 xa 2 (r]){-6(^' + HA) 2 + 2A^(2A - V) + 4A(5 + E')(ip' + HA)} 



// DAD , , 

d 4 xa 4 (7?) (m 2 s (A + AE — 3^) 2 - m 2 {3^ 2 + {AE) 2 - 2^AE + -=- + A 2 }) J .(86) 



B is non-dynamical, and for m\ ^ it is determined in terms of the other fields. For m\ = m 2 , 
A appears only linearly in the mass term. For the flat case H = and 0(77) = 1, this makes A a 
Lagrange multiplier and thus its variation gives rise to a constraint between the fields E and ip, 
leaving just one scalar propagating degree of freedom. In the de Sitter case, the result is the same, 
although this is not so obvious from the previous expression for the action until one substitutes 
the constraints. 

The variation with respect to A and B yields the constraints 

B (8T) 
a{r]) z mf 

= -2a(r 7 ) 2 m 2 (H(0 / - 3g) + A^) - a( V )*m 2 m 2 t (<l> - 3ij>) - SAW (RR) 

m 2 (m 2 s -m 2 )a(r]) 4 + 8An 2 -6m 2 a(r]) 2 n 2 ' 1 ' 
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where = AE. Let us first consider the kinetic part of the action, which after insertion of the 
constraints reads 

K = °^L{M l {r 1 W + M 2 ( V W 2 + MainW + ^)<A (89) 

where we have performed a partial integration to eliminate the term flip. The functions Mi(rj) 
are given by 

M M = 8A(m 2 -2m 2 s )a(v) 2 n 

u /; m 2 {m 2 - m 2 t )a(r)) 4 + 8AH 2 - 6m 2 a(r)) 2 H 2 ' 1 ; 

yr i \_ 2K ~ m\)a{rj)\\A - 3m 2 a(v) 2 ) , , 

2{V) m 2 (m 2 -m 2 )a(7]Y + 8An 2 - §m 2 a(r)) 2 H 2 ' 1 } 

M (ri) = 4m 2 {m 2 - m 2 )a{ V ) 4 

31 U m 2 (m 2 s - m 2 )a(r]) 4 + 8AH 2 - 6m 2 a{r]) 2 n 2 ' 1 ' 
-Am 2 a(r]) 2 H 2 

Ma{7] > = m 2 (m 2 - m 2 )a(r]y + 8AH 2 - 6m 2 t a{r}) 2 H 2 ' (93) 

A difference between the flat and the de Sitter backgrounds is that the coefficients M\ (rj) and 
M^{rj) cancel in the former case, and this automatically yields a kinetic term with a negative 
eigenvalue unless M 3 (rj) = which happens for the Fierz-Pauli combination m 2 = m 2 . The situ- 
ation in de Sitter is slightly more complicated. 



Let us now show that the previous kinetic term gives a positive contribution to the Hamiltonian 
in the range of parameters 

m 2 > 0, < m 2 - m 2 < 6H 2 . (94) 
Indeed, the kinetic term can written as 

'< - ^ + (mm - §£) r + mm (, + ^J) . m 

In the range (1911) . M^rj) and 4M 4 (r/)M 2 (^) — M 2 (rj) are positive. By Euler's theorem, the corre- 
sponding Hamiltonian 

R K = LV + IL^' - K, (96) 

is numerically equal to the two last terms in the Lagrangian, which are quadratic in generalized 
velocities, and hence it is positive definite. The second condition in (19~4"1) for a positive kinetic term 
reduces to the usual m 2 = m 2 for the Minkowski limit H = 0. For H > the endpoints of the 
interval are of different nature: the condition m 2 — m 2 > is necessary condition for positivity of 
M2 — MI/M4 at any value of the momentum, whereas the upper bound on the range of m 2 — m 2 
can be somewhat relaxed depending on the value of the momentum. Indeed, what we need is that 

m 2 — m 2 < 6H 2 



4A 

3a 2 m 2 1 



(97) 
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so the condition is considerably relaxed at wavelengths shorter than the inverse graviton mass. 



Once we have established the positivity of part of the Hamiltonian, let us see what happens 
to rest of it, namely to the potential part. This part will be given by 

V = K _ L = ^l!(M 5 (r / )0 2 + + M 7 (r^ 2 ), (98) 

where the coefficients are rather cumbersome and we omit them. Before proceeding, it should be 
noted that the Hamiltonian we are considering is time dependent, and hence not conserved. Its 
positivity and boundedness is a useful criterion only as long as we consider time-scales shorter 
than the expansion time, or energies larger than H. This is what we may call the adiabatic limit. 
Hence, let us assume that m s ,m t 3> H, even if their difference is much smaller m 2 — m 2 < H 2 , 
so that we can satisfy the positivity of the kinetic term as discussed above. We have checked that 
within this adiabatic limit, the potential V grows negative and unbounded below for — A/a 2 ^> m 2 . 
Instabilities at high momenta have been previously studied in [13], and they are just as bad as 
ghost instabilities. Unlike the case of tachyons, the phase space for instability is infinite and this 
yields infinite decay rates. 

If the masses m t and m s are small, of order of the expansion rate H, then we are outside of the 
adiabatic limit, and the Hamiltonian above is not a very useful indicator of stability. Instead, we 
should use a conserved charge associated to the time-like Killing vector for length scales smaller 
than the horizon [H]. Due to the existence of the cosmological scale, it is in principle possible 
(although by no means clear) that there may be some range 

m 2 (-Y >-A/a 2 >H 2 >m 2 , (99) 



(with a > 2), where this conserved charge is positive definite. The effective theory would then be 
well defined for momenta larger than H (corresponding to modes within the horizon), provided 
that the theory is cut-off at the energy scale m(H/m) a l 2 . We leave the study of this conserved 
charge for further research. We note, however, that we need a theory which is applicable to wave- 
lenghs much smaller than the horizon —A/a 2 ^> H 2 , where the adiabatic approximation should 
again be valid. We have checked that for —A/a 2 ^> H 2 > m 2 , the potential V grows negative 
and unbounded below, so the possibility of a range of the form (1991) where the conserved charge 
is positive does not look particularly promising. 

Finally, for the case m 2 = m 2 the analysis of the degrees of freedom has already been performed 
in another foliation in [15] (see also [46]). In our analysis for this case we find il^r?) = M^{rf) = 
and thus ip is not a propagating field. After varying the action with respect to ip we obtain a 
constraint which after substitution yields the Lagrangian 

/ 2 9mfr(77) 4 // H - 21m 2 a(ri) 2 fi H A + ±{p H - 6m 2 a(r/) 2 )A 2 + 2A 3 



(9m 2 t a(r]) 2 fXH - 6/i#A - 2A 2 ) 
where 2fi H = 2TL 2 — m 2 a(r]) 2 and 

3a(r]) 4 m 2 fiH 



T 



K + (9m 2 a(ri) 2 fi H - 6/i H A - 2A 2 
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This Lagrangian will be ghost-free and tachyon-free for fijj < 0. This reduces to the well known 
condition m 2 > 2H 2 [4"T] . 

Another case which differs from the usual Fierz-Pauli term and is still well defined is the 
case of Lorentz-breaking mass terms [T5], [16] . These mass terms arise from solutions of the form 
sketched in Section 2. A simple example would be given by two de Sitter metrics with A, = 
{/57 _1 , 7 _1 , 7" 1 , 7 -1 }. The interesting ranges propagating only one scalar or no scalars at all are 
the same as those in [T51 [T6], and for generic mass terms with two scalars, there is always a 
gradient instability at high energies as happens in the Lorentz preserving case [39J . 

5 Conclusions 

Bigravity is a natural arena where a non-linear theory of massive gravity can be formulated. This 
theory is invariant under the "diagonal" group of diffeomorphisms, under which both metrics 
transform. 

We have shown that if the interaction between the two metrics is non-derivative, we can always 
find exact solutions in the Schwarzschild-(A)dS family, with Lorentz breaking asymptotics. These 
are "Type I" solutions [2S] , in which both metrics cannot (generically) be brought to diagonal form 
in the same coordinate system. For a given interaction potential, the degree of Lorentz-breaking 
is an adjustable parameter. 

Perturbations around Lorentz-breaking bi-flat (or bi-de Sitter |39j) solutions lead to gravitons 
with Lorentz-breaking mass terms. Because of the invariance under diagonal diffeomorphisms 
[22] , mass terms with components are absent from the second order Lagrangian. This, in turn, 
leads to a well behaved theory of linearized perturbations [221 El], which is not afflicted by the 
vDVZ discontinuity. It is somewhat puzzling that in the linear theory, there are corrections to 
the Newtonian potential which are proportional to the square of the graviton mass and which 
grow linearly with the distance to the origin. On the other hand, as mentioned above, these 
theories admit the Schwarzschild metric as an exact solution for the same values of the parame- 
ters. Thus, the linearized solutions for a static spherically symmetric sources do not coincide with 
the linearization of the known vacuum solutions. This seems to indicate that this theory has a 
linearization instability such as the one which is found in other contexts [H], some of which are 
related to massive gravity and may have important phenomenological consequences [49] . Another 
possibility is that there may be other exact solutions which coincide with the linearized approxi- 
mation at large distances, and those may be the relevant ones which can be matched to spherically 
symmetric matter sources near the origin. This issue clearly deserves further investigation. 

We have also considered perturbations to solutions where both metrics are proportional to 
each other, focusing in the case of de Sitter. This has led us to consider generic mass terms for 
gravitons in de Sitter space, beyond the Fierz-Pauli case. We find that for the case of Lorentz- 
invariant mass terms, only the Fierz-Pauli combination is free from instabilities at high momenta 
—A/a 2 3> m? (and only for m 2 > 2H 2 ). For the Lorentz-breaking case, the situation is analogous 
to that in flat space [151 ESI [39] . 
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